Statistical analysis of a dynamical multifragmentation path 
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A microcanonical multifragmentation model (MMM) is used for investigating whether equilibra- 
tion really occurs in the dynamical evolution of two heavy ion collisions simulated via a stochastic 
mean field approach (SMF). The standard deviation function between the dynamically obtained 
freeze-out fragment distributions corresponding to the reaction ^^®Xe-|-^^®Sn at 32 MeV/u and the 
MMM ones corresponding to a wide range of mass, excitation energy, freeze-out volume and nuclear 
level density cut-off parameter shows a unique minimum. A distinct statistically equilibrated stage 
is identified in the dynamical evolution of the system. 
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It is known from more than 15 years that violent heavy 
ion collisions lead to an advanced disassembly of the com- 
pound system known under the name of nuclear fragmen- 
tation. Good agreements between various observables re- 
lated to the asymptotically resulted fragments and var- 
ious models assuming statistical equilibrium Q, |^ |^ Q 
lead to the conclusion that a huge part of the available 
phase space is populated during the fragmentation pro- 
cess. Hence the break-up of a statistically equilibrated 
nuclear source could be at the origin of fragment produc- 
tion. The source size, its excitation energy and its volume 
are thus quantities which can only be indirectly evaluated 
by comparisons between experimental data and statis- 
tical multifragmentation model predictions via a back- 
tracing procedure. However, the comparison process is 
complicated by the presence of several effects, such as 
pre-equilibrium particle emission, collective radial expan- 
sion (see e.g. 0), Coulomb propagation of the break-up 
primary fragments, secondary particle emissions. So it 
is difficult to ascertain whether these quantities really 
correspond to the source properties obtained at a given 
time during the fragmentation path. Moreover, there are 
intrinsic dynamical characteristics such as the freeze-out 
specific time, directly related to parameters such as the 
level density cut-off parameter, r (see Ref. contribut- 
ing to the weights of the system statistical ensemble. In 
the absence of any direct information about the freeze- 
out this last parameter has to be employed as a fitting 
parameter in a statistical model. In Rcf. [g, a very good 
agreement between statistical fragment distribution pre- 
dictions and experimental data was obtained assuming 
T = 9 MeV. However, one has to have direct access to 
the freeze-out events in order to unambiguously decide 
on the value of such parameters and, if an equilibrated 
source exists, to find its location in space and time. This 
task can be achieved using "freeze-out" information from 
a dynamical model. 

This letter investigates whether statistical equilibra- 
tion occurs in the dynamical path of two heavy ion colli- 
sion and, if so, which are the corresponding freczc-out pa- 
rameters. To this aim, the "freeze-out" data of a stochas- 
tic mean field (SMF) approach is analyzed via a sharp 



microcanonical multifragmentation model (MMM) 

We use the stochastic mean-field approach introduced 
in ref.Q. According to this theory, the fragmentation 
process is dominated by the growth of volume (spinodal) 
and surface instabilities encountered during the expan- 
sion phase of the considered excited systems. The dy- 
namical evolution of the system is still described in terms 
of the one-body distribution function (mean-field descrip- 
tion), however this function experiences a stochastic evo- 
lution, in response to the action of the fluctuation term. 
The amplitude of the stochastic term incorporated in the 
treatment is essentially determined by the degree of ther- 
mal agitation present in the system. Hence fluctuations 
provide the seeds of the formation of fragments, whose 
characteristics are related to the properties of the most 
unstable collective modes Q. 

In the model fluctuations are implemented only in 
r space. Within the assumption of local thermal equi- 
librium, the kinetic fluctuations typical of a Fermi gas 
are projected on density fluctuations. Then fluctuations 
are propagated by the unstable mean-fleld, leading to the 
disassembly of the system. 

The MMM model concerns the disassembly of a sta- 
tistically equifibrated nuclear source {A, Z, E, V) (i.e. 
the source is defined by the parameters: mass number, 
atomic number, excitation energy and freeze-out volume 
respectively). Its basic hypothesis is equal probability 
between all configurations C : {Ai, Zi, £i,ri,pi, i = 
1,...,N} (the mass number, the atomic number, the 
excitation energy, the position and the momentum of 
each fragment i of the configuration C, composed of N 
fragments) which are subject to standard microcanoni- 
cal constraints: J^i'^i = J2i^i = J2iPi = 0' 
X Pi = 0, -E - constant. The fragment level den- 
sity (entering the statistical weight of a configuration) 
is of Fermi-gas type adjusted with the cut-off factor 
exp(— e/r): /9(e) = /9o(e) exp(— e/r) 0. The above factor 
counts for the dramatic decrease of the lifetime of frag- 
ment excited states respective to the freeze-out specific 
time as the excitation energy increases (i.e. earlier freeze- 
outs should correspond to larger values of r). The model 
can work within two freeze-out hypotheses: (1) fragments 
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are treated as hard spheres placed into a spherical freeze- 
out recipient and are not allowed to overlap each-other or 
the recipient wall; (2) fragments may be deformed and a 
corresponding free- volume expression is approaching the 
integration over fragment positions 6]. Though more 
schematic, hypothesis (2) seems more adequate for the 
present study, so we use it herein. Further a Metropolis- 
type simulation is employed for determining the average 
value of any system observable (see Refs. U for more de- 
tails). While the model includes secondary deexcitation, 
this stage is not needed for this investigation so we will 
use here only its primary decay stage. 

Using MMM we investigate whether the "primary 
events" produced by the stochastic mean-field approach 
as a result of ^^^Xe-|--'^^^Sn at 32 MeV/u reaction may cor- 
respond to the statistical equilibration of the compound 
system. We consider only very central collisions of the 
above reaction. According to the dynamical simulations 
performed in Ref.0, it is observed that, after the initial 
collisional shock, the system expands towards low den- 
sities entering the unstable region of the nuclear matter 
phase diagram (after about 100 frn/c from the beginning 
of the reaction). Then fragments are formed through 
spinodal decomposition. The dynamical freeze-out time 
is defined as the time when the fragment formation pro- 
cess is over. Hence average fragment multiplicities and 
distributions do not evolve anymore. For the reaction 
considered this time is 240 fm/c. 

Our aim is to investigate whether fragment distribu- 
tions are compatible with the statistical phase space oc- 
cupancy, as predicted by MMM. For washing-up pre- 
equilibrium effects which should appear in the dynamical 
simulation, only intermediate mass fragments (IMF) (i.e. 
fragments with Z > 3) are selected. Therefore, all com- 
parisons between MMM and stochastic mean-field results 
are to be restricted to IMF's. Due to the large Coulomb 
repulsion among primary fragments it is reasonable to 
assume that the largest uncertainty in the equilibrated 
source estimated from the dynamical approach concerns 
the volume. In other words while real equilibration may 
occur at a different volume, than the one corresponding 
to the "dynamical" data, we assume that fragment sizes 
and excitation energies are roughly preserved. Therefore, 
we will fit the fragment size distributions and their inter- 
nal excitation energy but not the volume. The best fit 
can be found by minimizing the following error function: 

£ = {3[f {{Abound)) + fiiZbound))] 
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where (•) stands for average, Abound and Zbound are 
the bound mass and charge (sum of the mass number 
and, respectively, atomic number of all IMF's from a 
given event), Njmf is the number of IMF's, ejMF is 
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FIG. 1: Contour plots of the error function £ [see eq. (1)]: in 
the (V/Vb, E) plane corresponding to A = 210 (upper panel); 
in the {A,V/Vo) plane corresponding to E — 5.7 MeV/u 
(lower panel). Darker regions correspond to smaller £; units 
are relative. 

the fragment excitation energy per nucleon and Zmaxi 
with i = 1,2,3 are the largest, second largest and third 
largest charge from one fragmentation event. Further, 
f{x) = \2{xs —Xd)/{xs + Xd)\, where the indexes s and 
d stand for "statistic" and "dynamic" and is the ab- 
solute value of X. One may observe in eq. (1) that 
the bound mass and charge terms are overweighted with 
an (arbitrary) factor 3. Of course, this does not in- 
fluence the significance of the error function, £, which 
will vanish for a perfect agreement between the s and 
d results. On the other side, this will put stronger con- 
straints on source size identification. For finding the best 
agreement between statistical and dynamical results we 
variate the MMM parameters A, E, V and r in suffi- 
ciently wide ranges thus constructing a four-dimensional 
grid. The ranges arc A : [195,230], E : [3,7.5] MeV/u, 
V/Vq : [1.5, 9.5], T = 12, 16, oo MeV. The source is consid- 
ered to have the A/Z ratio of the ^^^Xe-|-^^^Sn reaction. 
Considerable computing power was involved for obtain- 
ing this result in reasonable time. An absolute minimum 
of the error function, £, was found at A = 210, Z = 87, 
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TABLE I: Evolution of Emin (see text) with r. 
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V/Vo = 3.4, E = 5.7 MeV/u, r = oo. Cuts in £ corre- 
sponding to A — 210 and E = 5.7 MeV/u are represented 
in the upper part and respectively the lower part of Fig. 
1. The evolution of fmin (i-e. the minimum £ corre- 
sponding to a given t) with r is represented in Tabled 
From the above mentioned figure one can see that the 
global minimum is clearly determined and there are no 
secondary minima. 

The statistical source corresponding to the minimum 
value off yields the following results: (Abound) = 199.03, 
(Zbound) = 83.74, (eiMp) = 4.21 MeV/u, (Z„,axi) = 
42.34, (Z„,ax2) = 24.35, (Z^naxs) = 11.4. These are to 
be compared with the corresponding dynamical results: 
{Abound) = 199, (Zbound) = 84, (eiMp) = 4.3 MeV/u, 

(^maxl) = 41.95, (Z,„ax2) = 22.5, (Z,„ax3) = 13.3. NotC 

the excellent agreement for all considered observables, 
proving the very good quality of the fit. That means 
that the fragment size related features obtained dynam- 
ically have been reproduced by the primary decay stage 
of the MMM model. This can be seen in Fig. 3, as well, 
where the dynamical Z and Njmf distributions fit per- 
fectly the MMM ones. Due to the fact that the fragment 
size distributions and fragment excitations were used for 
constructing £ these results are natural, given the small 
value of the obtained minimum. At this stage the ques- 
tion still remains: even if fragment size distributions and 
excitation energies of the fragments are very well repro- 
duced, do the dynamically formed fragments come from 
an equilibrated source with freeze-out volume V — 3.4Vb, 
as predicted by MMM? 

In the present work the freeze-out volume V is the vol- 
ume of the smallest sphere which totally includes all frag- 
ments. We denote by V the volume of the smallest sphere 
which totally includes all fragments and has the center 
located in the center of mass of the system. Obviously, 

V > V. The "dynamical" events have (^{V /Vo)imf'^ = 

9.08; the statistical ones have ^(T^/Vb)/Af_F^ = 4.93 (the 
IMF index indicates that we refer to the volume occu- 
pied by IMF fragments). This means that equilibration, 
as predicted by MMM, may have occured at an earlier 
time compared to the considered dynamical events, i.e. 
pre-fragments already appear at an earlier time and they 
are actually equilibrated inside a smaller volume. 

One can simply test this hypothesis: one just has to 
propagate the fragments in their mutual Coulomb field 
from the freeze-out positioning as generated by MMM 
up to {V /Vo)iMF = 9.08 (i.e. the value corresponding to 
the "dynamical" events) and then compare "dynamical" 
and "statistical" fragment kinetic energies and positions. 
However, in performing such a comparison, one has to 
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FIG. 2: Dynamical charge (upper panel) and number of IMF 
distributions (lower panel) in comparisons with the statistical 
ones. Statistical results are represented by open squares; 
histogram (upper panel) and stars (lower panel) corespond 
to the dynamical ones. 

guess the initial fragment velocities, that could have a 
collective component. Flow is easily accontable for in 
MMM 1^]. The best reproduction of the dynamical re- 
sults is obtained for a flow energy equal to zero. (Note 
that this result is in agreement with the dynamical simu- 
lations where the initial expansion energy present at early 
times is reduced by the strong interaction between pre- 
fragments and finally flow becomes close to zero.) The 
time of the Coulomb propagation is around 100 fm/c. 
The comparison is presented in Fig. 4: the fragment 
average kinetic energies versus charge is represented in 
the upper panel; the radial distribution of the fragment 
with largest, second largest and third largest charge is 
represented in the lower panel. An excellent agreement 
is observed between "dynamical" and "statistical" data 
for both observables indicating the physical consistency 
of the scenario used in the present statistical analysis. 
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FIG. 3: Upper panel: "statistical" (open squares) and 
"dynamical" (histogram) fragment average kinetic energy 
versus charge. Lower panel: "statistical" (solid lines) and 
"dynamical" (dashed lines) radial distribution of fragments 
with largest (peak "1"), second largest (peak "2") and third 
largest (peak "3") charge in one fragmentation event. The 

plot corresponds to (^{V /Vo)imf^ ~ 9.08. 

It means that the dynamically generated events may re- 
ally originate from a statistically equilibrated source with 
V = 3.4Vb. 

Some considerations are now in order. All fragment 
properties at the dynamical freeze-out appear to be al- 
most univocally linked to a statistical emission from a 
more compact equilibrated source at an earlier time- 
step, equilibrium freeze-out^ that eventually expands on 
Coulomb trajectories. This time (i ~ 140 — 150/m/c) 



can be interpreted as when the dynamical evolution of 
the system reaches a stage of global equilibrium. 

The equilibrium freeze-out time range appears just in- 
termediate between the beginning of the spinodal decom- 
position and the final fragment configuration. At that 
time the leading unstable modes are well established and 
some pre-fragments, in strong interaction with the rest of 
the system, can be recognized. In this situation a statisti- 
cal fragmentation picture appears adequate, apart some 
minor effects, as the few "equal fragment size" events, 
relics of the space structure of the primordial spinodal 
instability [9] . However, due to the expansion of the sys- 
tem, the dynamical freeze-out configuration, where frag- 
ments are well formed and separated, corresponds to a 
larger volume with respect to this equilibrium freeze- out. 

Since in dynamical simulations fragment formation is 
coupled with monopole oscillations of the source, so the 
system may recontract in some events or evolve towards 
more dilute configurations in some other cases it is 
obvious that we have to deal with the fluctuation of the 
freeze-out volume from event to event. On the other 
hand, in MMM calculations while V/V[) is fixed, V /Vq is 
fluctuating. Whether or not the fluctuations of V /Vq are 
sufficient to mimic the "dynamical" ones or one should 
explicitly include fluctuations of V/Vq it remains an open 
question. Since larger fluctuations are reflected into 
Coulomb effects it would be interesting for future studies 
to compare the kinetic energy distributions around the 
average of Fig. 4 (top) for each charge. 

From the nuclear thermodynamics point of view an 
important conclusion is that one deals with small freeze- 
out volumes (3.4 Vb) and highly excited primary frag- 
ments (around 4.3 MeV/u), i.e. hot fragmentation at 
small freezeout volumes. As it is well known, volume is a 
key variable for locating the system in the phase diagram 
of the nuclear matter and a lot of work aiming to identify 
this quantity is in progress |llj| . 

Summarizing, along the fragmentation path, as de- 
scribed by a dynamical model, a huge part of the avail- 
able phase space is filled. Statistical equilibration occurs 
in a pre- fragment stage of the system. To our knowledge, 
this is the first time when in dynamical paths of violent 
heavy ion collisions a stage of statistical equilibrium of 
the compound system is proved. 
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